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Starting from microscopic models of space-time foam, based on brane universes propagating in 
bulk space-times populated by DO- brane defects ( "D-particles" ) , we arrive at effective actions used 
by a low-energy observer on the brane world to describe his/her observations of the Universe. These 
actions include, apart from the metric tensor field, also scalar (dilaton) and vector fields, the latter 
describing the interactions of low-energy matter on the brane world with the recoiling point-like 
space-time defect (D-particle). The vector field is proportional to the recoil velocity of the D- 
particle and as such it satisfies a certain constraint. The vector breaks locally Lorentz invariance, 
which however is assumed to be conserved on average in a space-time foam situation, involving 
the interaction of matter with populations of D-particle defects. In this paper we clarify the role 
of fluctuations of the vector field on structure formation and galactic growth. In particular we 
demonstrate that, already at the end of the radiation era, the (constrained) vector field associated 
with the recoil of the defects provides the seeds for a growing mode in the evolution of the Universe. 
Such a growing mode survives during the matter dominated era, provided the variance of the D- 
particle recoil velocities on the brane is larger than a critical value. We note that in this model, as a 
result of specific properties of D-brane dynamics in the bulk, there is no issue of overdosing the brane 
Universe for large defect densities. Thus, in these models, the presence of defects may be associated 
with large-structure formation. Although our string inspired models do have (conventional, from a 
particle physics point of view) dark matter components, nevertheless it is interesting that the role 
of "extra" dark matter is also provided by the population of massive defects. This is consistent 
with the weakly interacting character of the D-particle defects, which predominantly interact only 
gravitationally. 



I. INTRODUCTION AND MOTIVATION 



Relativistic modified gravity theories {e.g. models in Refs. [UH]) have been presented originally as field theoretic 
alternatives to dark matter models, offering support to the heuristic Modified Newtonian Dynamics approach [3 . 
Although, at least in their simplest form, such models may be ruled out by means of recent precision astrophysical 
measurements using gravitational lensing 0], nevertheless some of their features may characterise microscopic theories 
of quantum gravity, in harmonic co-existence with substantial dark matter components, thereby avoiding the above- 
mentioned stringent constraints. 

For instance, it has been shown in Ref. [5], that the bi-metric nature of the TeVeS (Tensor- Vector-Scalar) model, 
as well as the existence of a Lorentz-violating vector field A^, do characterise the low-energy limit of certain string- 
theory models of space-time foam [6l [7] , involving three-brane ( "D3-brane" , a three-space-dimensional Dirichlet brane) 
worlds [8j [9] embedded in higher-dimensional (bulk) space-times punctured by DO-brane ("D-particle") defects. In 
such models, the D3-brane representing our Universe, which may be obtained from appropriate compactification of 
a higher-dimensional brane, moves in the bulk and in this way the D-particles cross it, resembling — from a D3- 
brane observer viewpoint — flashing "on" and "off" foamy structures ("D-foam"). There are topologically non-trivial 
interactions of open string states, attached on the D3-brane and representing ordinary matter, with such defects, 
involving splitting of the initial string and creation of intermediate string states, stretched between the D-particle and 
the D3-brane universe, exhibiting length oscillations. Such processes result in local distortions of the neighbouring 
space-time [6] 

Ak 

g^u=r]^u+h„ v , h i=Vi=g B j^, (1) 

to leading order in the recoil velocities «j<l, due to recoil of the D-particle defect to conserve energy and momentum 
during the scattering. In the above formula, g s is the (weak) string coupling, M s is the string mass scale, which is in 



2 



general different from the four-dimensional Planck scale Mp = 1.2. x 10 19 GeV, and Vi is the recoil three-velocity of 
the D-particle, with Afc^ the momentum transfer of the open string state. 

Due to electric charge conservation, only electrically neutral matter interactions with the D-particles 1 are allowed. 
In this sense, there is a naturally induced bi-metric structure in the model, given that it is only the electrically neutral 
excitations of the low-energy effective theory that predominantly feel the induced metric Eq. (JlJ. Moreover, like in 
TeVeS models, there exist vector field structures, provided by the recoil velocities of the D-particles, which break 
Lorentz invariance locally. 

As pointed out by the authors of Ref. [TT] , it is the vector field of the TeVeS theory that they were working with, 
which plays the crucial role in reproducing the observed power spectrum by generating a cosmological instability that 
produces the large cosmic structures of today's Universe. It is important to note that the arguments on the role of 
the vector field in reproducing large scale structures and the correct phenomenology seem to be generic, in the sense 
that they do not depend on the detailed Lagrangian of the original TeVeS theory proposed by Bekenstein [Jj . More 
precisely, it is only the basic important features of the theory, namely the existence of the two metrics (bi-metric 
theory) and of the (Lorentz-violating) "aether-like" vector field, that appear to be important in this respect. 

The basic aim of this paper is to extend the link between the above mentioned features and the deduced cosmological 
perturbation theory, in the framework of some specific backgrounds of string theory proposed some time ago [5] as 
consistent candidates for a quantum space-time foam background in string theory. Specifically, we shall construct a 
low-energy effective action, describing (partly though) the low-energy aspects of such a stringy space-time foam, and 
argue that it acquires the form of a modified gravity theory, with scalar and vector fields arising naturally, as a result 
of the interaction of matter with the D-particle space-time defects in the foam. Then, we shall consider the equations 
for cosmological perturbations, solve them to linear order for small perturbations, and argue that the vector fields in 
the theory play a non- negligible role in structure formation, in the spirit of Ref. |llj . However, we stress that such 
theories do not provide alternative to dark matter, given that the latter exists naturally in the superstring-inspired 
models we are working with, represented by stable superpartners to the standard model matter that characterises 
the low-energy excitations of the models. In this sense, our approach is entirely different from models which are 
alternatives to dark matter [TJ [2J. Nevertheless, the D-particle defects themselves play a role analogous to dark 
matter, given that their recoil furnishes the Universe with an effective vector field, whose perturbations in late eras 
of the Universe (radiation and matter dominated) do posses growing modes and can participate non-trivially in large 
structure formation. Thus the massive D-particles can, in addition to baryonic and other type of matter, contribute 
to galaxy and galactic cluster formation. 

To set-up the framework we shall use, we briefly review in Section [TT] the brane D-foam model that will serve as 
our microscopic framework for discussing modified gravity in the low energy limit. First we present the underlying 
formalism of the world-sheet deformation that describes the interaction of an open string state with a recoiling D- 
particle and derive the associated metric distortion of the neighbouring space-time, due to the recoil of the defect. We 
explain in detail the emergence of a vector "gauge" field as a result of this interaction, which will play a crucial 
role in our analysis in providing us with the growing mode that can lead to structure formation. The vector fields 
are proportional to the recoil velocities of the D-particles, and as such they satisfy a certain type of constraint (cf. 
Eq. ( 52 1), which is crucial for the appearance of the growing mode. We discuss in Section III more general background 
configurations for the dilaton, graviton and gauge fields, in a way consistent with the conformal invariance of the world- 
sheet of the string. We first present the Dirac-Born-Infeld (DBI) action on 3-brane worlds and we then proceed with 
space-time curvature correction terms in the DBI action. In Sect ion [TV] we study the four-dimensional induced effective 
action on a D3-brane world for the gravitational, dilaton and gauge fields. This is the action that we subsequently 
use to analyse the role of vector perturbations in structure formation. In Section [V] we derive the equations of motion 
for the vector, metric and dilaton fields and consider their cosmological perturbations. In Section VI we discuss the 
appearance of a growing mode in both the radiation and matter dominated eras of the brane Universe, provided a 
sufficiently large variance of recoil velocities of the D-particles exists. In this sense, the defects act as a dark matter 
component, which is in addition to the conventional (from a particle physics viewpoint) dark matter components of 
string effective models, e.g. supersymmetric partners of matter or graviton fields. Finally in Section |VII| we give our 
conclusions and outlook. 

We use the following conventions: for the metric wc use the signature (—,+,+,+) and the Ricmann curvature 
tensor is defined as i2g 7f = d s T a pi +T x ^T a xs - ( 7 < ► S). 



This is the case of type IIA string theory, which allows point-like branes. In the phenomenologically more realistic type IIB string 
models, there are no DO-brane configurations allowed. In such a case, one can still construct consistent D-foam models |10l in which the 
role of the D-particle is played by D3-branes wrapped up around small three cycles. In such a case, electrically charged excitations do 
interact non-trivially with the D-foam, but these interactions are significantly suppressed compared with those of the neutral excitations. 
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FIG. 1: Schematic representation of a generic D-particle space-time foam model. The model of Ref. [B], which acts as a 
prototype of a D-foam, involves two stacks of D8-branes, each stack being attached to an orientifold plane. 



II. REVIEW OF D-FOAM INSPIRED MODIFIED GRAVITY EFFECTIVE THEORIES 



The prototype model that we shall use as a microscopic motivation of our considerations in this work is the so- 
called D-brane model of space-time "foam" , in which our Universe is represented as a (compactified) three-brane world 
moving in a higher-dimensional bulk space-time which is populated by point-like DO-branc defects (D-particles) . Such 
a model is appropriate for type IIA string theory and the associated supergravity effective field theory in the bulk 
and was first presented in Ref. [BJ. There is a relative motion of the defects with respect to the branes, which implies 
that, as the brane universe moves in the bulk, the D-particles cross it, so that from the point of view of an observer on 
the brane they look like "flashing on and off" foamy structures (hence the terminology "D-particle foam" model). In 
string perturbation theory, it can be shown that for an adiabatic relative motion between a D-particle defect and the 
brane, with the defect moving perpendicularly to the brane, the induced potential on the brane, due to open strings 
stretched between the D-particle and the brane, depends on both the relative velocity v± and the distance of the bulk 
defect from the brane: 



(long \2 

Vno S -m = + A ^ L ,r»Va> (2) 



However, a D-particle close to the D3-brane (compactified D8), at a distance r' <C \[a! , moving adiabatically in the 
perpendicular direction with a velocity Uj^ ort ,will induce the potential : 



II short \ 2 

-i;short _ na \ V ± I fn\ 
y D0-D8 — 10 ,3 W 



This difference in sign, implies that, one can arrange for the densities of far away and nearby bulk D-particles, which 
are not in general homogeneous, to be such that the total contribution to the brane world's vacuum energy is always 
sub-critical, so that issues such as over- closure of the Universe by a significant population of D-particle defects can 
be avoided [7]. Consequently, the astrophysical data can only constrain the total energy density on the brane due to 
the defects, that is the algebraic sum of the above two contributions in Eqs. ( pj), pj ). Hence, the density of defects on 
the brane remains unconstrained. We shall come back to this point in Section jVI| 

These considerations imply that, at least at early eras of the Universe, we can plausibly assume that the brane 
Universe passes through bulk regions which are densely populated by D-particlcs, in such a way that there is a 
sufficiently high density of D-particles on the brane, such that the particle excitations, represented by open strings 
with their ends attached on the brane, propagate in a "merfwm" of such D-particle defects. Topologically non-trivial 
interactions between the defects and matter strings can occur, involving the capture and splitting of open strings 
by the defects and creation of stretched strings between the defects and the brane during such processes [BJ [7] (c/. 



4 



Fig. [T]) . The capture process is represented from a world-sheet point of view as an impulse which leads to a metric 
deformation of the neighbouring space-time due to the recoil of the D-particle. The induced metric distortions depend 
on both the space-time coordinates and the momentum transfer of the matter particle during its scattering with 
the defect. The recoil of the defect breaks locally Lorentz invariance of the space-time as a result of its momentum 
direction, and in the effective low-energy limit is described as a vector field A^, which should not be confused with 
the electromagnetic field. On average, over large populations of defects Lorentz invariance may be restored in the 
sense of having a zero vacuum expectation value of the vector field, but non-zero variances. The purpose of this 
work is to examine the effective low-energy target space action describing the interaction of matter excitations with 
the D-particles, through the field and discuss in this context the potential role of the perturbations of the vector 
field in inducing a growing mode that could participate in large-scale structure formation at redshifts of order z ~ 1. 
As a result of charge conservation, charged open string excitations cannot interact with the (neutral) defects via 
the above-mentioned splitting. Only electrically neutral excitations are therefore interacting predominantly with the 
D-particle foam. Cosmologically it is therefore the neutrinos and photons that feel mostly the effects of the foam at 
these early eras, and hence it is these particles that can find themselves propagating in the background of recoiling 
DO-branes, and hence of the associated vector field A^. We shall show here that D-particles can play the role of a 
dark matter component, in addition to either supersymmetric partners of standard model excitations or neutrinos. 

We commence our discussion towards a formulation of an effective low-energy action for this type of stringy inter- 
actions, including the dynamics of the vector field itself, by first reviewing briefly the underlying formalism of the 
world-sheet deformation that describes the (impulse) interaction of an open string state with a recoiling D-particle [H] . 
The pertinent vertex operator on the world-sheet boundary <9£, within the world-sheet cr-model reads: 

v;Zf c = / u t x°e e (x°)d n x* , (4) 

JdT. 

where X° is the target-space time coordinate, obeying Neumann boundary conditions, and X 1 (i = f,2,3) are the 
spatial coordinates on the brane, obeying Dirichlet boundary conditions. The notation d n denotes normal derivatives 
on the world-sheet. By Ui we denote the recoil relativistic three- velocity of the D-particle {v,i = jVi = jg s Aki/M s , 
with 7 = (1 — v 2 )^ 1 / 2 the corresponding Lorentz factor, and Aki the momentum transfer of the stringy-matter 
excitation). The operator 

<3 e (X°) = -i r J^- e ^ x ° with e -> 0+ , (5) 



is a regularised operator. It can be shown [12] that for e -2 = 2\n(A/a 2 ), where A/a 2 is the area of the world-sheet in 
units of a world-sheet Ultra- Violet (UV) cut-off length scale a, the vertex operator in Eq. Q satisfies an appropriate 
logarithmic conformal field theory. The operator has an anomalous dimension — e 2 /2 — > + in the Infra-Red (IR) 
fixed point A /a 2 —> 00, where we shall work from now on [5]. 

The impulse operator in Eq. Q is written in the co-moving frame of the recoil D-particle. In the frame of an 
observer on the D3-brane, the covariantised form of this expression should be used instead [5]. We can write this as 
a total world-sheet derivative in the following way: 

/ d a {u l ,u u X L '@ e {u p X p )d a XL l ) with e^0+, (6) 

where is a four-velocity vector satisfying: 

= -1 , (7) 

and a = 1,2 is a world-sheet index. Note that above, using the Neumann boundary conditions, d n X° — 0, for the 
temporal coordinate A , thereby implying an ambiguity in Eq. Q when written as a total world-sheet derivative, 
allowed us to write the covariant form Eq. ([6]). 

As discussed in detail in Ref. [5], where we refer the interested reader for details, the bulk operator Eq. ^ is 
conformal on the world-sheet in the sense of having conformal dimension two. However, the target-space metric 
deformation implied by the vertex Eq. ([6| does not connect smoothly with the flat Minkowski space-time before the 
capture of the open string by the D-particle. Indeed, by expanding the derivative in Eq. ^ and ignoring terms of 
the form d a d a X fl , that can be eliminated on-shell in the world-sheet theory, we recognise a contribution similar to 
the free sigma-model action propagating in a curved background g^, namely J" s g fJlU d a X fJ, d a X' / . In this sense, there 
is a target-space deformation due to the D-particle recoil: 

Sg^ = u^UvQ^+lupX 11 ). (8) 
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To ensure a smooth connection with Minkowski space-time before the collision/capture, one can couple to this de- 
formed a- model a linear dilaton background term of the form [5]: 

d 2 i(f>R (2) , 4> = u^X" , (9) 

where <i 2 £ is the covariant measure over the curved world-sheet of curvature R^ 2 \ Coupling the vertex operator 
Eq. ^ to factors does not affect the conformal invariance, since the conformal dimension of in the presence of 
the linear-dilaton background Eq. ^ is u^(u^ — u^) — 0. Hence the complete vertex operator for recoil, compatible 
with conformal invariance at the IR fixed-point e — > + , turns out to be [5]: 

V rccoil = f e^d a (u^u u X"Q^ a+ (u p XP)d a X^ . (10) 



Acting with the derivative in the above formula, we can see that the deformed target-space metric has the form 

5g„ v = e*« M u„e e _ >0 + (upX"). (11) 
However, to ensure a smooth connection with the flat metric at the origin of the boosted time, u p X p = 0, the above 



metric deformation Eq. (11 1 must also contain extra dilaton terms, such that the induced metric is of the form: 



neutral matter _ _ ( <j> _ p -<t>\ 

= t^X" . (12) 

Note that the correction corresponds to an operator on the world-sheet of the string with zero conformal dimension, 
— Ufj,(— — u^) + 2 = 2(w Al w Al + 1) = , due to Eq. This latter deformation leads to departure from criticality of 
the associated cr-model and thus Liouville dressing is required [13] . 

To this end [5], we first notice that the linear dilaton implies a sub-critical string with Q 2 = u^u^ = — 1 < 0. 
This can become conformal if one uses a space-like Liouville mode |13| p to "dress" the above-mentioned metric 
deformations by multiplication with exponential operators e QiP (with i = 1,2); a, denote the Liouville "anomalous" 
dimensions. In the presence of the world-sheet background charges in the (p,X^) extended target space time of the 
Liouville-dressed world-sheet theory, induce world-sheet curvature terms of the form 



in the cr-model action. The conformal dimension of these operators is 0^(0^ + 1) for each i = 1,2. Restoration of 
conformal invariance requires that the total conformal dimension of the Liouville-dressed deformations is (1,1) in the 
(holomorphic, anti-holomorphic) world-sheet sectors. It is straightforward then to observe that the following dressed 
operators are conformal, amounting to the choices on = ±|Q| = ±1 in the respective Liouville anomalous dimensions, 

V&*»d)«A«* = (e-^-P - e u » x * + ?) dX x dX»u x u v . (13) 

These imply a dressed target-space-time metric in the extended space-time (p, X) of the form: 

^neutral matter, dressed^ X ) = ^ + L,Hp,X) _ e -*(/>,A^ u ^ > 

9 Pl i = o , 
9pp = +1 ( 

${p,X) = -u^-p. (14) 

The extra space-like Liouville mode may thus be given the physical interpretation of a bulk spatial dimension, in 
which recoil of the defects does not take place. Here, our brane space-time is located at, say, p = 0. Note that this 
is just one example of a consistent conformal theory. In general, one may consider non-trivial cr-model metrics G pv , 
in which case the associated dilatons will have a more complicated space-time dependence. In this respect, one can 
discuss the effects of D-foam recoil in realistic Friedman-Lemaitre-Robertson- Walker (FLRW) backgrounds, which 
constitutes the subject of the present paper. 

With the above considerations in mind, we take from now on the following form of the deformed target-space metric: 

^neutral matter = Q ^ + {p4 > _ with ^ v = Q, . . . 4 , 

4> = $(0,^) with /i,i/ = Q,...4, (15) 
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due to D-particle recoil during their interaction with neutral matter. Thus, we assume a three (spatial)- dimensional 
(D3-)brane universe from now on for brevity (the latter may be obtained by appropriate compactification of higher 
dimensional branes), ignoring bulk physics as far as the recoil of D-particles is concerned, which we assume to be 
dominant only on the brane world. 

Finally, we remark that the recoil velocity field is elevated to a dynamical one in this approach, as being part 
of a gauge background field in the D-particle recoil, which obeys non-trivial dynamics [5 . This can be seen in 
two equivalent ways. In the first one, we use the fact that, upon T-duality (which exchanges Neumann and Dirichlet 
boundary conditions and is assumed to be an exact symmetry of the underlying string theory), the (non-covariantised) 
vertex operator of recoil Eq. Q is related to the background gauge potential deformation [T^] : 

M = Ul X Q <d(X°) , (16) 

assuming |12) . without loss of generality, a time axial gauge A = 0. This background corresponds for X° > (after 
the impulse) to a constant "electric" -type field Ei = F oi = u i: where F^ = d^A u is the Maxwell field strength of the 
gauge field. This vector field describes the average effects of the interaction of neutral matter with the background of 
the recoiling D-particle and should not be confused with the electromagnetic field of the standard model; it is a new 
degree of freedom associated with the back reaction of the D-matter onto space time. 

Alternatively, without invoking T-duality, one may use the world-sheet version of Stokes' theorem, to write down 
the boundary recoil deformation Eq. Q (before coupling to the dilaton/Liouville) as a bulk deformation [14] : 



■^impulse 



27ra' 
1 

47ra' 



7 / d^{2u i )e a ^X° 



e £ (x°) + x°s £ (x 



d a X i , (17) 



where e a p is the world-sheet Levi-Civita antisymmetric symbol and 8 e {X ) is an e-regularised 5-function. For relatively 
large times after the impact at X° — (which we assume here for our phenomenological purposes), this is equivalent 
to a deformation describing an open string propagating in an antisymmetric .B^-background corresponding to an 
external constant in target-space "electric" field, 

B 0i = F oi ~ m and By = for X° > , (18) 

where the X°S(X°) terms in the argument of the electric field yield vanishing contributions in the large time limit, 
and hence are ignored from now on. In this approach, we note that the presence of the B-field leads to mixed-type 
boundary conditions for open strings on the boundary <9£ of world-sheet surfaces with the topology of a disc; we 
consider here for concreteness and relevance to our discussion below, that: 

9^d n X v + B^d T X v \ d v = , (19) 



where g M „ denotes the metric of target space-time. Formally, the operator Eq. (17) is conformal on the world-sheet, in 
the limit e — > + 2 in flat space-times. The generalisation to curved space, andm particular to FLRW backgrounds, 
will be discussed in subsequent sections. 

This conformal gauge field background is assumed to exist together with the metric and dilaton backgrounds, 
discussed previously. The dynamics of the gauge field describes physics in the so-called open string channel, whilst 
the dilaton and metric backgrounds describe the effects of the recoiling D-particles in the closed string channel. The 
above considerations summarise therefore the features of the space-time foam model which motivate the inclusion of 
vector and scalar fields together with gravitational tensor fields. 



The reader should notice a formal similarity of the metric Eq. ( 15 ) with the corresponding one in TeVeS models [T], 
including, in addition to the graviton tensor field g^, also scalar (dilaton cf>) and vector (recoil velocity-related gauge 
Ap) fields. However, in our case, the interpretation is entirely different from the TeVeS models, and moreover, 
given that we are dealing with ordinary (from a phenomenological point of view) supersymmetric strings, there are 
natural candidates for dark matter in our models (the lightest (stable) supersymmetric partners to the standard model 
excitations), co-existing with the scalar and vector structures. Thus, our models are not meant to provide alternative 



2 Notice that to express the gauge deformation in the frame of a D3-observer, it is not sufficient to replace X° — > u^X^ , in the argument 
of the O e _ > Q+(X°), but also to transform the recoil three-velocity vector Ui. However, for our purposes this is not necessary. 
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to dark matter scenarios. Nevertheless the existence of dilaton and vector fields can affect certain cosmological features 
of the string universe, which shall be explored in this article. 

In a cosmological context, as already mentioned, one needs to reconsider the conformal invariance conditions for 
metrics of FLRW type that depend on the cosmic time. These conditions are equivalent to equations of motion for 
the various target-space fields that are derived from a low-energy string effective action. We next proceed therefore 
to discuss the effective low-energy action, which can describe in general quantum fluctuations of the tensor, scalar 
and vector backgrounds and in subsequent sections we shall discuss background solutions to the respective equations 
of motion and perturbations around them. The background solution we shall deal with in this work is characterised 
by a constant dilation but non trivial FLRW metrics and recoil vector fields adapted to this geometry, which will 
generalise the flat space solution ( 16 1 mentioned above. 



III. LOW-ENERGY EFFECTIVE LAGRANGIAN FOR D-FOAM MODELS: A PROPOSAL 



Here, wc shall generalise the arguments of the previous section to discuss more general background configurations 
for the dilaton, graviton and gauge fields, in a way consistent with the conformal invariance of the world-sheet of the 
string. We would like to search for cosmological backgrounds to be applied in our cosmological tests of the model. 
In particular, we will investigate the role of vector fields to galaxy growth. We shall construct first an effective low- 
energy field theory for a graviton, dilaton and vector field background in the context of the model described above 
and find configurations for the background fields that are consistent with the conformal invariance conditions of the 
cr-model theory, generalising the analysis of Ref. [5]. Then, we shall analyse the resulting modifications of gravity 
in a cosmological context, in particular we will study their effects in galactic growth and structure formation. The 
requirement of conformal invariance of the string theory on the world-sheet is equivalent to the requirement of the 
background fields satisfying equations of motion arising from the target space-time effective action. As the vector-field 
part of the action used to describe open strings ending on D-branes is well-known and can be expressed in a closed 
form to all orders in the string Regge slope a' (DBI action) |T5j, we are able to start from it and then study its 
modifications to include the description of gravity in the bulk and on the brane. 



A. Dirac-Born-Infeld action on the 3-brane worlds 



We commence our analysis from the gauge field effective Lagrangian in four-dimensional flat space-time, representing 
the case of propagating open strings in gauge-field backgrounds, with their ends attached on a D3-brane. The gauge 
field backgrounds and the open strings are not propagating in the bulk. In the open string channel, it is known that 
such a Lagrangian can be obtained to all orders in the a' expansion, and has the DBI form [15] : 

£dbi = --T 3 Jdct 4 (r^ + 2na'F^) , (20) 

where T3 is the three-dimensional brane tension and F^ V (A) is the (Abelian) field strength. It can be readily shown 
that a constant electric field strength background, as is the case of the recoil-velocity gauge field for times after the 
impulse, X° > 0, is a consistent solution of the gauge field equations of motion derived from this action, and thus a 
consistent (world-sheet conformal invariant) configuration for the string. 

In the presence of non-trivial dilatons, the string coupling g s becomes g s oe^ (where <? s o = e << ^ > is a space-time 
constant). Moreover, in non-Minkowski backgrounds, the tensor rj^ v should be replaced by the background metric 



g^ v and thus the DBI Lagrangian Eq. ( 20 1 acquires the form: 



£dbi - -— T 3 e-*Jdet 4 (g^ + 2i:a'F^) + ... , (21) 

where the . . . include other terms, which are non-trivial in the case where the field strength is non-constant, and 
also include curvature terms that give the graviton field dynamics. In the closed string channel, of course, the latter 
terms are just the Einstein-scalar curvature terms to lowest order in derivatives for the gravitational sector. However, 
an interesting question arises as to whether one may discover curvature terms coupled directly to the DBI action, as 
a result of induced metrics in the open string channel. It is the point of this section to attempt to address such a 
question. 

First of all we remark that the four-dimensional DBI action can be expanded in derivatives, as appropriate for a 
low-energy approximation, compared to the string scale M s = 1/ sfa 1 , as follows |15) : 

det 4 (ff M „ + 2na'F^) = det 4 .g [l + {2^a') 2 h - (2W) 4 /f ] , (22) 



where 



h = -g^g up F^F Xp , I 2 
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Upon replacement in the Lagrangian Eq. (21 1, one obtains: 

V dct 4.9 



C 



DBI 



+ l (2na') 2 e-^h - J(2W) 4 e^- if + • • • + £ 



fJsQ 



.9s() 



(23) 



(24) 



where the . . . indicate higher derivative terms, and Cq indicates closed-string sector gravity contributions. Since 
gravity is allowed to propagate in the bulk, the four-dimensional form of the gravitational terms can be obtained from 
a projection onto the D3-brane of such terms. In the cr-model frame the latter are of the form: 



1/(6) 

9%M! 



(25) 



where is a volume element of extra (compactified) dimensions (in units of n/o'), A2 is an effective cosmological 
constant term, which may arise from compactification, and tf> is a four-dimensional dilaton kinetic term. The signature 
of A2 depends on the details of the theory. We define the four-dimensional bulk-induced gravitational constant kq &s: 



1 



y(e) 



(26) 



Notice that in this model there are two contributions to the dilaton-induced dark energy: one from the open string 
sector, as a result of the 



9s0 



terms in Eq. ( 24 1 , and another one from 



y(e) 



e- 2 *A 2 = l(<t>) 



terms. Combining Eqs. (24 1, (251, the overall dark energy contributions would then come from 



A 



total dark 



(27) 



(28) 



(29) 



Since in our approach, Ai is manifestly of fixed sign (positive, for positive tension D3-branes T3 > 0), the overall sign 
of the induced four-dimensional dark energy A tota i dark energy depends on the relative strength of its constituents. On 
noting the different scaling of the two types of vacuum energy contributions Ai, A with the string coupling g s oe^, we 
observe that it is possible to have a positive dark energy on the four-dimensional world, for positive Ai(0) + A(</>), 
with the Ai contribution being sub-dominant for weak string couplings at late times. This is the assumption we shall 
make when analysing the consequences of the Eq. ( 24 ) in cosmic structure (galaxy) formation. For other phases of the 



brane universe, the sign of A((f>) may change. In fact, in the microscopic models of D-foam [6] [7] we are considering 
here, this sign depends on the ratio of the densities of nearby-to-far-away bulk D-particle defects, with respect to the 
brane universe. Indeed, the attractive flux forces between brane and bulk defects moving in directions transverse to 
the brane world, have specific dependence on the relative distance, which are such that for nearby (far away) bulk 
D-particles (lying at distances shorter (longer) than the string length from the brane universe) there are negative 
(positive) energy contributions to the brane vacuum potential energy. We shall come back to this important point 
later on in our article. 



The effective actions in Eqs. (24 1 and (25) do not give the complete action that couples dynamical gravity to the 



recoil vector field; some further modifications need to be done to the Dirac-Born-Infeld part, which we shall now 
discuss. 



B. Space-time curvature correction terms in the Dirac-Born-Infeld action 



We first notice that, in addition to the closed-string sector gravitational field, the effective action of D-foam may 
contain induced target-space curvature contributions, that couple to the open-string sector, in particular to the DBI 
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Lagrangian Eq. (21). Such a situation is known to characterise simple examples of Dp-branes (p stands for the 
dimensionality of the D-brane), fluctuating in the bulk space [16]. In fact, the authors of Ref. [16] represented the 
essential features of such fluctuating branes by considering toy examples of world-sheet actions for a Dp-brane with 
tachyonic-type deformations of the form 



1 

Ana' 



dX M dX 



C 

Ana' 



or. 



(30) 



where X M with M = 0, ...,D- 1, are the string coordinates. In this notation, the index [i takes values from to 
p and corresponds to the p + 1 coordinates that obey Neumann boundary conditions, while the rest obey Dirichlet 
boundary conditions; the latter are denoted by X 1 . The limit c — > oo, which corresponds to a conformally invariant 
point of the theory, ensures that the end-points of the strings are confined to move on the hyper-surfacc defined by 
Y l . In this limit, a space-time effective action for the Y l s can be identified with the partition functional of c and Y % : 



S(Y l ) = lim Z{c,Y l ) = ( dX l dX^e 

c— >oo J 



The authors of Ref. |16) split the string fields in the usual way in a zero mode plus perturbations, 



X 



M 



X 



M 



(31) 



(32) 



and Taylor expanded Y t (X) around Y % = Y l (x). They then integrated out of the partition function the contribution of 
X, and, after expanding to orders in (x l — Y l )K % ^ u: with K % = d^d^Y 1 being the extrinsic curvature, they integrated 
the quadratic and quartic terms and arrived at 



Z(c,Y) = Z(c) / dxVdetG 



l-a'((2)R(G)+(D(a' 2 ) 



(33) 



where C(2) = tt 2 /6 and G„ v = ^ + d^Y l d v Y\ 

In the flat-brane limit of Ref. |16j . the induced metric G is the only gravitational contribution. One, however, may 
generalise the background G to an arbitrary metric, g, which should then be a consistent solution to the classical 
equations of motion for the gravitational field, obtained from this modified gravity action. The presence of the 
curvature corrections in the partition function Eq. ( 33 ) of the open string, imply a target-space effective action of the 
form 



5' 



9s 



J d^x^J— det 



<J 



a'^-R + 0(a' 2 ) 
6 v ' 



(34) 



which is just the ordinary Einstein-gravity theory (to lowest order in derivatives), in the presence of a positive 
cosmological constant T^,/g s . From this point of view, the gravitational (Planck) constant for this theory is k — 
g s /(a'Tz). However, the reader should also notice that there are curvature contributions to the effective action, 



coming from the closed string sector, of the form Eq. (25). In the case of constant dilatons (as is the case of the 
model of Ref. [161 ), this will lead to additional contributions to the effective gravitational constant in four-dimensions, 
coming from the closed string sector, while in the case of time-dependent dilatons, of interest to us here, one would 
obtain field-varying gravitational couplings. We shall come back to this point in the next section. 

For the moment we note that coupling a U(l) gauge field background to the model, Eq. (30), corresponds to adding 
the following boundary term to the action S^: 



V2 



h L A ^ xv ^ 2vh L - 



(35) 



where we have used Stokes's theorem, and the indices a, j3 denote world-sheet indices. The \j2na' comes from 
dimensional considerations. For constant field strength backgrounds this is formally equivalent to an antisymmetric 
5^-tensor background. In general, this is not the case. This addition changes the effective action Eq. (33) in the 
two following ways. 

Firstly, the determinant of the (induced) metric changes as 



y/Actg \f&ct(g + 2na'F) = y/det(h) , 



(36) 
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where we consider a general metric background g, according to our previous discussion, and not only the induced 
metric G. As we have discussed in the previous sub-section, the determinant given in Eq. (36) is the one that appears 

is defined as 



in the DBI action in a curved background, Eq. (21 ). In fact, the tensor h^ v 



2-kolF.. 



(37) 



and its inverse is denoted by W" v , with W" v h vp = S^ p , Secondly, as it is well-known from matching with string- 
amplitude calculations |17j . in addition to the curvature terms in Eq. (34), there are contributions involving the 
gauge field and its target-space derivatives (in the general case of non-constant field-strength backgrounds). In the 
superstring theories we are interested in (despite the fact that we are only considering the bosonic parts of the 
pertinent target-space effective actions), it is known that such corrections start from structures with four-derivatives 
acting on the field strengths, i.e. of the form 18J 



5i 



DBI 



3s 
9s 



d 4 a; v /det(/i) 



1 — a'—R + C urv(a' 2 ) + O ((ddF, OFdF) 2 ) 



(38) 



where O curv (a / ) indicates terms involving curvature squared terms. All such higher-derivative terms will be ignored 
in our low-energy considerations. 



IV. THE FOUR-DIMENSIONAL INDUCED EFFECTIVE ACTION 



Following the above discussion, we therefore use the four-dimensional effective action (on a D3-brane world) 



Jeff 4dim 



d 4 x 



— e -0 yj- det (g + 2na'F) [ 1 - aR{g) 
9s0 \ 



~ ge -^ A A + - ^ger 2 '*' R(g) + 0((d^) 2 ) 



(39) 



(where g — det(g)) for the gravitational field g^ u , the dilaton field 4> an d the gauge field A^, with field strength 



F^u — d^Ajy 



dyA^ on a D3-brane. Note that k q is defined in Eq. (26), and a — a'C,(2) — a'ir 2 /6 in the example of 



Ref. |16j , however we wish to keep the discussion general and so we use from now on the positive constant a > as a 
parameter of our model; it has dimensions of length squared. In this respect, we also assumed a non-constant dilaton, 
thereby writing the string coupling as g s — g s o^ , with g s o < 1 a (perturbatively) weak constant string coupling. The 
terms 0((d(f>) 2 ) denote kinetic terms of the dilaton, which shall not play an important role in our analysis in this 
work, for reasons that will be clarified below; hence we do not need to write them down explicitly. 



Using Eq. ( 22 ) , taking into account that for the background solutions configurations pertinent to the effects of the 
D-particle foam only "electric-field" type components F$i are non-trivial (i.e., 1% = in Eq. (|23|) and expanding the 



square root in Eq. (39), while keeping up to four-derivative-order terms, we finally arrive at the following form of the 



four-dimensional effective action on the D3-brane world 



~cff 4dim 



d 4 x 



9s0 



K 



^2to' e^F^F^ + a- F ^2Tra' e'^F^ F uv R(g) 



9s0 



3 - 

a — e 



1 



e" 2 ^ R(g) + O ((def,) 2 



(40) 



In view of the mis- match between open and closed string sectors, we do observe from Eq. (40) that the coefficients of 
the Einstein scalar curvature term R are not of the usual Brans-Dicke type. 



We shall base our subsequent analysis upon Eq. ( 40 1 . We may proceed by redefining the space-time metric in such 
a way that the bulk-induced Einstein curvature term (1/kq) e -2 ^ R(g) acquires a canonically normalised term without 
dilaton couplings, that is we redefine: 



9^f ^ 9{iv e ^9^ 



(41) 



which defines the so-called Einstein- frame metric |19j . In our conventions, such a redefinition implies for our four- 
dimensional brane world: 



a\ - e 2 * 



R(g)- 



6 



: d„ (V=ff 9^ M + 6 g^d^ d v <f> 



(42) 
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Upon combining Eq. (41 1 with a canonical normalisation of the kinetic terms of the gauge fields, by absorbing the 
factor {2-Ka' e~^°T^/ g s oy/ 2 into a redefinition of the gauge potential A^, we finally obtain the following effective action 
in the Einstein frame S^ e 4dim (for brevity from now on we omit the tilde notation from the Einstein metric g^v): 



ScS 4dim — ' • 



"4 



-F^F^e-^R{g) + (a**- + -) R(g) - § a e 2 V F^)(d^) + O' ( W) 2 ) 

4 V 9so K oJ 1 



S m , (43) 



up to terms proportional to the square of dilaton gradients 9 M </>, denoted by O' ((d<j>) 2 ), which we shall ignore, as we 
have already mentioned. The reader should have noticed, however, that terms containing one derivative of the dilaton 
have been kept, since they will yield non-trivial contributions to the dilaton equation of motion even in the case of 



constant dilaton, we are considering. In Eq. (43) we have added the matter sector action S m which we do not specify. 
Notice that for effectively constant dilaton <f> = <j>Q we are interested in here, the gravitational constant (that is the 
coefficient of the Einstein term) reads 

- = — +a— e* . (44) 
K k g s0 

This is the action we shall use from now on in order to analyse the effects of the vector perturbations in structure 
formation in this stringy universe. Notice that, for the sake of brevity, we used above the same symbol F for the field 
strength of the canonically normalised gauge fields. 

Before proceeding to discuss the role of perturbations, let us discuss first the background configurations about 
which we shall perturb. First we may assume that at early epochs of the universe, which will not be the subject 
of our study here, there is a de Sitter phase of the brane universe, in which the population of D-particle defects is 
suffi cien tly dilute, so that any recoil vector field strength F contributions to the effective action are negligible. From 



Eq. (43) we observe that apart from the term — (Ts/g s o) e 3 ^, the rest of the terms in the effective action have the 
form of a low-energy gravitational action coming from the closed sector of string theories, within graviton, dilaton 
(<fi) and tachyon (T(x)) backgrounds, the latter one corresponding to initial cosmological instabilities. Such de Sitter 
phases have been argued to be consistent with world-sheet conformal invariance to all orders in a' Ref. [2D]. In fact, 
the de Sitter phase corresponds to a relatively brief cosmological era, during which the amplitudes of the (classical) 
tachyon and dilaton fields are proportional to each other, T(x) + <f)(x) = 0. Once such a field alignment is lost, due 
to a decay in time of the tachyon instability, the universe exits dynamically the de Sitter phase. It is important to 
notice that in the solution of Ref. [2DJ , the dilaton in the Einstein frame also assumes in such a phase a logarithmic 
Robertson- Walker-frame time dependence, <fi = fiahit, <po < 0. This implies that the afore- me ntioned open-string 



term — (T^/gso) e 3 ^ is negligible compared to the rest of the terms in the effective action Eq. (43) and thus it can be 
approximated by the one derived from closed strings. This justifies the existence of an approximately de Sitter phase 
in our context. 

At late eras, for instance those corresponding to red-shifts around z ~ 2, relevant for structure formation which we 
are interested in, the D-particle populations, especially near galactic centres, are significant, thereby leading to the 
presence of recoil vector field ^-contributions to the effective action. In such eras, the background configuration for the 
D-particle recoil operators which satisfies the (logarithmic) conformal algebra on the world-sheet of the string [2T] , and 
is therefore a consistent solution of the equations of motion of the target-space effective action, involves approximately 
constant dilatons and spatially flat Robertson- Walker metric backgrounds in the Einstein frame. 



This, in conjunction with Eq. (15), gives us that 



4> = (f> Q = const , g 00 = -1 + 2sinh((/> )u 2 l , g iS = a 2 (t)(5y + 2sinh(0 o )'U i u i ) . (45) 

However, in what follows we shall redefine the coordinates in such a way so that the standard FLRW metric is used 
as our space-time background in order to make contact with observations. 

The configuration for the D-particle recoil operator in such backgrounds reads 3 : 

V= f g^{t)y 3 {t)Q{t - t imp&ct )d n x l , (46) 

JdY, 



3 For completeness we note at this stage that, even if we worked with the metric Eq. | |45| , the factors depending on sinh(</>o)fi^ can 
be ignored below, because they give rise to terms cubic in Ui or higher. However, by performing the above-mentioned coordinate 
transformations, we make sure that no extra factors involving the dilaton 0o would hamper the comparison with the data. 
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where d n denotes the normal derivative on the world-sheet boundary 9S, t is the cosmic time, gij{t) = a 2 (t) is the 
spatial part of the FLRW metric and 



1 - 2p 



^-impact j with t ^> ^impact -> 

(47) 



is the D-particle recoil trajectory, with v the D-particle recoil three velocity in the co-moving frame and ti mpa ct the 
moment of impact of the matter string with a single D-particle. Let us note that this expression stems from the 
geodesic trajectories of a recoiling D-particle in the metric background |21) . 

In our study we consider the interaction of strings with populations of D-particles, in which case ii mpac t is an averaged 
time, considered in what follows as a free constant. Let us, for convenience, absorb the constant impact /(I — 2p) in 
the definition of velocities, which also absorbs the normalisation factors (27ra'e~^°T3/g B o) 1 ^ 2 employed above in order 
to normalise the Maxwell terms in Eq. (43). This is equivalent to considering a vector gauge field background with 
spatial components of the form: 

Mt) = -a 2 {t) 5 ijV j = -a(t)S ljV l phys) , (48) 

where w ( p h ys ) = a (t) v ^ > 3 — 1)2,3 is the physical ("local") recoil velocity. The above analysis holds for low three 
velocities; otherwise, one should replace in Eq. (48) the three velocity by the spatial components of the four- vector of 



the velocity, u l = jv 1 with 7 = 1/ Wl — v t vigij. 

The background Eq. (48) corresponds to a field strength of "electric field" type: 

F 0i (t) = -2a{t)a{t) Ul . (49) 

where the over-dot denotes the derivative with respect to cosmic time t. 

We remind the reader that in this formalism [12 , conformal invariance of the stringy cr-model does not restrict 
the form of the temporal component of the four- vector field A^, apart from the fact that homogeneity and isotropy 
require it to depend only on the cosmic time, Ao(t). This freedom allows us to covariance the background vector field 
as follows: 

A» = -a 2 {t)u ll = -a{t)uf* s . (50) 



In the Roberson- Walker background of Eq. (45), the physical ("local") four velocity u^ ys = aifiju^ obeys the 
Minkowski- flat constraint: 

u (phys) u (phys) _ _| constant | _ _Eil 2 7rQ! / e-^ < , (51) 

gsO 

from which it follows that the vector field has the following time-like constraint in our Robertson- Walker background: 

AuA^ = -^27ra'e-*° < . (52) 
9s0 

Unlike the case of Ref. pQ, however, in our model the vector field has non-trivial space-lik e co mponents, proportional 
to a 2 {i)ui. When we consider perturbations, we shall maintain the above constraint, Eq. (52), which is implemented 
in a path integral via an appropriate Lagrange multiplier A(x) term in the effective action, namely 

SW-Lagrangc = J d^X^X^) ^# + M 2to V*J , (53) 

which notably does not couple to the dilaton field. The above should be added to the right hand side of Eq. (43). 

We next remark that, upon averaging over foam populations, one assumes the following Gaussian stochastically 
fluctuating configurations: 

< Ui >= and <C Ujltj >= a 2 Sij , (54) 

where <C • • • S> indicates appropriate averages over D-particles populations in the foam. Thus, we observe that only 
upon averaging over the isotropic foam, we obtain vanishing of the spatial components of the A^ field and of the 
"electric field", Eq. (49), namely 



< Ai >= and < F Qi >= 



(55) 
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Nevertheless, in our case caution should be exercised to take these averages only at the very end of the computations, 
since even powers of Ui yield non-zero population averages. The parameter ctq is free in our model, albeit small, to 
be constrained by the data. In general, Ctq is cosmic time dependent since it is a function of the (bulk space) density 
of the foam, which may vary for different cosmological eras [Hl[7]. In general, assuming a no- force condition among 
the D-particles [7J, which characterises D-foam models derived from super- membrane theories [BJ, implies that the 
density of foam scales as "dust" with the cosmological scale factor, a 3 (t), so that to a good approximation over the 
time period of structure formation we are interested in here, we may take 



al(t) = L^J ( p = constant > . (56) 



Here the factor 2-Ka'e~^° \T^\/g S Q comes from the constants which were previously absorbed into the definitions of the 
recoil velocities. 

Some important remarks are due at this juncture. The above considerations are in agreement with the isotropy 
and homogeneity of the observed Universe. However, for completeness we should mention that when considering the 
effects of our model at galactic scales, it is possible that certain inhomogeneities in the density of the D-particles may 
occur in such a way that D-particle populations dominate the haloes of galaxies, thus playing a role analogous to dark 
matter. Indeed, due to their localised nature (as contrasted with the space filling D-brane universes ), D-particles 
may also be considered as massive excitations of the vacuum with mass of order M s /g s . Since there is a no-force 
condition among them, such particles are ideal to play the role of super- weak dark matter (termed D-matter |22j). 
and if they are sufficiently light, as is the case of low string scales M s of order TeV, and g s = 0(1), they may have 
implications at colliders, where they can be produced by the collision of Standard Model particles [3 [22]. However, 
as discussed in Ref. [7J, and reviewed briefly above, in our model of D-brane universes propagating in bulk spaces 
punctured by D-particle populations, there are contributions to the vacuum energy on our brane world by the bulk 
D-particles, of mixed signature, that depend on the relative distance of the D-particle from the brane world. Negative 
(positive) energy contributions are due to bulk D-particles that lie at distances smaller (larger) than the string scale 
from the brane. As a consequence of such contributions, the density of D-particle populations that are trapped on 
our brane world is not restricted by the requirement of avoiding over-closure of the Universe and can thus be an 
arbitrary phenomenological parameter. In this sense, scenarios in which various mass density profiles for D-matter 
can exist in the haloes of galaxies, playing a role analogous to dark matter, are compatible with current cosmology. 



When we adapt our model above, Eq. (43 1, to such a situation, we observe that, in view of our coupling of the vector 
recoil field to the curvature of space-time, averaging over populations of D-particles in the galactic haloes, will result, 
at a background level, in extra contributions of the recoil term <C e _3c ^ F^ v F^ v ^> to an "effective" gravitational 
"constant" , proportional to erg, which in inhomogeneous scenarios may depend on the particular galaxy: 

' a « e-^ Ffj, v F» v » , (57) 



«eff 8 7rG off K 4 



where k is defined in Eq. (44) and plays the role of the four-dimensional Newton's constant Gn in our stringy frame- 
work. Such an effective gravitational constant will lead to Einstein equations coupling the matter and gravitational 
systems with a space and time varying effective G e s which can then affect cosmological considerations. Although such 
scenarios are worth pursuing further, especially as providing alternatives to conventional dark matter (in view of the 
absence of any concrete particle physics evidence for dark matter at present), nevertheless in the current discussion 
we shall restrict ourselves to homogeneous and isotropic densities of D-particles on our brane universe, which satisfy 



Eq. ( 56 ) . Our aim in this paper is to investigate the effect of the recoil vector field to the growth of galaxies and not 
to provide realistic alternatives to dark matter scenarios. Our super-membrane effective field theory has its own dark 
matter, provided, for instance, by the supersymmctric partners of the low-energy effective action. Nevertheless, we 
hope to come back to such interesting variants of our model in a future work. 



In the next section we proceed to consider perturbations around the background solution, Eq. (45), and discuss 
their relevance to structure growth. To compare our analysis with that of Ref. [llj , on the importance of the vector 
field for structure growth, we neglect in what follows fluctuations of the dilaton field, assuming them to be suppressed 
as compared to the gravitational and vector perturbations. This is consistent with the phenomenology of low-energy 
physics at late eras of the universe, given that such dilaton fluctuations would affect the values of the string coupling 
e^, and through it of all the coupling constants of the string- inspired effective theory, which would not be acceptable 
from a particle physics viewpoint. 
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V. BACKGROUND CONFIGURATION AND EQUATIONS OF MOTION 



By varying the effective action given in Eq. (43) with respect to the vector, metric and dilaton fields we obtain 
respectively 4 : 



ae 



R)]'' u + 2X(x)A ll =0 



(58) 



4 no g s o 



1 



\{x){A a A a + JZM 2 We-*°) 



9s0 



9jw 
2 



9s0 



-20o 



g aX F, lX F va {l-ae- 2 ^R) 



g^ 2 [F a0 F a p] - V M V„ [F^F a0 ] \=T£- X(x)A^ A v 



(59) 



and 



;e 2,p0 \7 2 (F afs F 



a0\ 



A 
! — t 

K 



,200 



— e 00 

ftrO 



1 



FavF^" e 



-20o 



i?( 3 ) = 



(60) 



where «o is defined in Eq. (|26j), is the matter stress tensor and we recall that a = o/£(2) = a'ir 2 /6 > in the 
D-brane example of Ref. |loj. Above, we have assumed for simplicity that the dilaton coupling to the matter action 
is not dominant (or equivalently that the dilaton couples only through exponential couplings of the form exp(ji<p), 
with 7j constants, to the various matter fields, such that upon considering the corresponding dilaton variations, one 
obtains terms proportional to the Lagrangians of the various mattcr-spccics i that vanish on shell). As our primary 
purpose here is to investigate the role of the vector perturbations on the structure formation and growth, we believe 
that such an assumption is reasonable. 



Before actually solving the system of equations ( 58 ) , ( 59 ) and ( 60 ) for the galactic growth and structure formation 



epoch of the universe, we consider it as instructive to discuss the emergence of a de Sitter phase in this system, which 
as we shall argue may indeed characterise late eras in the evolution of this stringy universe, thereby making our model 
consistent with the current cosmological observations indicating the dawn of an accelerating de Sitter evolution era 
once again in the history of our observable Universe. 



A. Condensates of the recoil vector field and late-epoch phases of the D-foam universe 



To demonstrate the existence of a de Sitter era we first remark that the dynamics of the system of the recoil 



vector field is described by a Dirac-Born-Infeld Lagrangian given in Eq. (21 ). For such lagrangians, the higher order 
interaction terms among the vector field strength may lead to condensates, in analogy with the gluon condensates 
of Quantum Chromodynamics. In fact, such an assumption has been made in Ref. |23| . for generic Abelian flux 
fields that characterise D-brane excitations, with Dirac-Born-Infeld world volume actions. The analysis in that work 
indicated then that under certain plausible assumptions on the dominance of the quantum effects (over classical, 
thermodynamical) on the formation of the condensates, it is possible to obtain an equation of state for the vector 
Dirac-Born-Infeld system that resembles that of a de Sitter phase, i.e. a cosmological constant type with w ~ — 1. 
Let us first review briefly, for instructive purposes, those results. Consider a generic Dirac-Born-Infeld field, with 



Lagrangian Eq. (20) on a D3-brane world volume. In general, there are two contributions to the vector field strength 



F^ v condensates and one may come from purely quantum vacuum effects, 

(F^ FHvac = 5(t) , (F^ F*Hvac - )§(*) , (61) 



4 We set to zero the terms involving derivatives of the dilaton only at the end of the variations, since the latter is assumed constant <f>o 
and absorbing the constant factors e~^ ^ 2 into a normalisation of the gauge potential A^ , so that in this case the kinetic term of the 
recoil gauge field assumes a canonically normalised Maxwell form. The constraint Eq. 152b , is then modified appropriately, however the 
dilaton equation remains unaffected, since the dilaton does not couple to the constraint term, as already mentioned. 
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where F* indicates the dual tensor, and the condensates may in general depend on time, but are constant on spatial 
hyper-surfaces. The isometry structure of the spatial hyper-surfaces lead the authors of Ref. .23] to further assume 
the following, which we also adopt here: 



(FouF ») 



a t (t) 



-9oo 



(FioF®) 



Os(t) 



9ij 



(62) 



where i,j spatial indices on the three-dimensional volume of the D3-brane and we have the relations a t + 3a s = Aa. 
In addition to the quantum effects, one also has classical thermodynamical effects on the energy density and pressure 
of the Dirac-Born-Infeld fluid, which are obtained by averaging over the spatial volume. If one decomposes the field 



strength into "electric-field", Foi = Ei, and "magnetic-field" Bi 
classical effects contributions on the condensates 123 : 



£ijk F 3 components, one may specify these 



(F „ F ") = (J2 El) , (Ft o F]) = - (E { Ej) + 2(B t B s ) 



(63) 



On making the further (natural) assumption [23 (Ei Ej) = (Bi Bj) — C gij/3, we observe that the total contribution 
of both classical and quantum effects to the vacuum condensates can then be expressed by the relations: 



a t — a — AC , a s — a + - C 



(64) 



Computing the stress tensor of the Dirac-Born-Infeld fluid, Eq. (20), for the case of the total condensates one arrives 
at the following expressions for energy density p DBI and pressure p DBl 23 : 




p DBI = 




A = 



(65) 



Notice that in our case (Eq. (20)) we did not subtract the determinant of the metric 
action 



/—g in our definition of the 
This explains the difference 
, by constant terms of 



as this plays the role of a cosmological constant contribution, which we keep 
between the right-hand-side of the above expressions and the corresponding ones in Ref. 
the form A/2 and —A/2 for p° B1 and p DBI , respectively. 

The quantum condensate a and /3 have not been specified. The only restrictions come from the positivity of the 
corresponding quantities inside the square root in the Dirac-Born-Infeld action, e.g. in the denominators of Eq. (65), 
which imply the following condition between the various quantum condensates: 



> -1 
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(66) 



From Eq. (64), it becomes immediately apparent that on the one hand, if quantum effects are the dominant ones, with 
a > C, then at — a s ~ a, and hence from Eq. (65) we obtain an equation of state of cosmological constant (de Sitter 
vacuum) type, namely u> DBI ~ — 1 independently of the exact form of the quantum condensate (assuming of course it 
exists, which is a question that probably cannot be addressed in a generic manner, as it requires specific properties 
of the brane action). On the other hand, as demonstrated in Ref. [23], when classical effects dominate, with C ^> a, 
then Eq. 



([64]) implies a t — — 3a s ~ —AC and thus Eq. (65) leads to an ordinary relativistic fluid with positive energy 
density and pressure, p DBI ~ (A/2)C > and p DBI = (l73)p DBI ~ (A/2)C > 0. 

Hence, the effects of the classical (thermodynamical) contributions to the condensates are to cause the equation of 
state of the Dirac-Born-Infeld fluid to deviate from the value u>dbi — —1 by an amount that depends on the relative 
strength of those effect with respect to the quantum ones. 

In the context of our study here, the Dirac-Born-Infeld vector field has a microscopic origin, as describing the 
dynamics of the recoil degrees of freedom of the D-particles under their interaction with electrically neutral matter. 
In contrast to the case of Ref. US] , here we have the additional couplings of the space-time curvature to the Dirac-Born- 



Infeld action, Eq. (21), whose effects have been ignored in Ref. [23]. Nevertheless, as we shall argue, for sufficiently 



small values of the condensate, one can still argue on a late de Sitter era, provided again the quantum vacuum effects 
dominate the formation of the condensate. 

At this point we would like to emphasise two important remarks. Firstly, in our case the classical effects on the 
condensates are just the stochastic effects obtained upon averaging over statistical populations of D-particles, Eqs. ( 54 ) , 
(|55| and ( 56 1 . However, the quantum effects have a genuine quantum nature due to the quantum-gravitational 



fluctuations of the D-particle defects in the absence of any matter. The latter are also described by the stretching of 
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virtual strings between the D-particle and the brane world but they are non-perturbative, and can only be computed in 
the strongly coupled regime of the target space action of the D-particle, therefore only within the context of M-theory 
and this is not yet known. Nevertheless, such effects do exist and in what follows we shall adopt a phenomenological 
approach and simply assume (as in the case of Ref. |23| ) that they can condense, when the D-particle populations are 
sufficiently dense on the brane world, and can lead to the quantum Dirac-Born-Infeld type condensates, Eqs. (61 1, 
p). 



Secondly, in contrast to the generic arguments of Ref. [53] , in the case we are studying here, there is the constraint 
Eq. (52), implemented in an effective action treatment through the Lagrange multiplier term, Eq. (53). There are 



two distinct phases in the problem, as a result of the quantum treatment of the constraint; we discuss them below. 

(a) The phase in which the constraint is irrelevant, i.e. (A(x)) = 0, in which case the vector Dirac-Born-Infeld field 
is massless. This is the case for an era in which sufficiently dense populations of D-particles cross our brane universe 
and they condense due to quantum effects, which are dominant over the stochastic population effects, Eq. (54). In 



this case, the quantum effects on the recoil vector field condensate dominate and one may write, according to our 
previous discussion and taking into account the isometry structure Eq. ( 62 ) , that 



Vac 



a{t) 



9vv 



(67) 



where the quantum condensate a(t) can be even negative as we have seen previously (c/. Eq. (66 )). In this phase, the 



recoil vector field represents quantum fluctuations of the D-particles in the absence of any matter string excitations. 
The nomenclature "recoil" may therefore be misleading but we keep it for uniformity purposes. 

(b) The phase in which (X(x)) = Ai 2 /2 > 0, in which case the vector field becomes massive. This is the phase 
where contributions to the growth of galaxies will come from, and is assumed to characterise the era of structure 
formation. Quantum effects on possible condensates are assumed sub-dominant. One has the classical effects on the 



condensates Eqs. (49), (54), (55) and (56) that dominate in this case 



1. de Sitter phase and quantum vacuum condensates 



We concentrate on the phase (a), which as we shall demonstrate, is compatible with a late-era de Sitter phase for 
our Universe, even after including the coupling of the target space-time curvature to the Dirac-Born-Infeld terms] 
Eq. (39), which was ignored in the analysis of Ref. 23 . To this end, we concentrate on the case where the condensate 



a is constant in time, small in magnitude, while the condensate j3 — 0, so that an expansion up to order a in the 
effective action is sufficient. To this order our perturbative equations (58), (59) and (|60|) suffice. Consider the dilaton 



equation (60) first. Upon assuming the formation of quantum condensates, Eq. (|67[> , we observe that this equation 



T e 50o 



implies that the scalar curvature of the space-time reads 

R = -. ir- 

a V^f e 3 ^ - a) 

This can be a positive constant, as appropriate for a de Sitter space-time, provided 

2T 3 



6 

9sQ 



4A 

K 



a > 



7mi 



(68) 



(69) 



If we assume a de Sitter maximally symmetric form for the Riemann curvature tensor, corresponding to a Hubble 
constant, H\ = const > 0, namely, 



we obtain from Eq. ( 68 1 the following relation 



Hi 



R 
12 



1 



( IS 

V 3=0 



T 3 c 50o 



A 

3k 



> 



(70) 



From the graviton equation (59), ignoring the matter contributions TJJl, using Eq. (67), and considering the phase 
(A) =0we obtain 
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From Eqs. ( 70 ) and ( 71 ) we finally obtain the value of the condensate that guarantees a de Sitter geometry in this 



phase, namely 



4 g20o 



1 



2k. e <Po i A 

g _ ffsO «o 

T 3 P0O + -A_ 

2g s0 3 kq 



Thus, the condensate is proportional to the factor 



a n Q 7r 2 ( 5s0 e-^) 5 



(72) 



(73) 



using Eq. ( 26 ) and adopting the value of a that appears in the example of Ref . [TB] . This factor depends on the size 
of the compactification volume. 

For mathematical consistency of our solution, the condensate a and the curvature, i.e. the Hubble constant H\ 



(701 have to be sufficiently small, which is easily guaranteed from Eq. (72) for sufficiently large negative values of 



cpo- Such large negative values of the dilaton may characterise late epochs of the universe. For instance, in the linear 
(run-away) dilaton scenario, in terms of the FLRW time in the Einstein frame, the dilaton assumes large negative 
values, <j) ~ — lni 3> —1 for large times t — > oo, and its variation is suppressed by 4> ~ — V* 0' so ^ ma y be 
considered approximately constant. This could be a phase where the above conditions for the formation of small 
condensates are satisfied, in which case the universe will enter a late de Sitter phase. 

From our microscopic D-particle foam point of view, such a late phase may occur when matter in the brane is 
sufficiently dilute due to the brane world cosmic expansion, and the brane universe passes again a bulk area with 
sufficiently dense D-particle populations. The classical recoil effects of the D-particles, due to their interaction with 
matter strings on the brane, become then sub-dominant to their quantum vacuum fluctuation effects. In fact, as 
already mentioned previously, in our microscopic models of D-foam Refs. @[7], in such densely populated bulk areas, 
the bulk cosmological constant A may be negative, 



A = — | A | < for dense populations of nearby bulk D — particles 



(74) 



since the negative contributions to the brane vacuum energy from the nearby bulk D-particles dominate. It is not 
inconceivable then, although admittedly this requires some fine tuning, that the bulk density populations are such 
that in order of magnitude 



Q > -1 



e-g 



9s0 



2\A\ 
3k 



in which case Eq. (72) implies: 



a 



2e 



•24>u 



To 
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(75) 



(76) 



and the condensate a -C 1 (as required for consistency of the approach) for |A|/«o "C 1 and finite values of <j>Q. 

In such a case, the formation of the condensate may be understood as stabilising the brane vacuum. Indeed, as 
follows from Eq. (59), upon the formation of constant quantum vacuum condensates a, the dark energy terms in the 



brane effective action assume the form 



Brane Dark Energy ~ a + e 



2<?n 



— 
9s0 



(77) 



In the absence of a condensate, the brane vacuum energy would be negative, of order — |A|//to and his would indicate an 
instability of the va cuu m. The true stable vacuum of the theory would then be the one in w hich the condensate forms, 



with the value Eq. ( 76 ) , which implies a positive (de Sitter- type) vacuum energy, Eq. i77\, of order e 2 * 1 ^ > 0, for 
any value of <po. This would guarantee the validity of our arguments on a D-vacuum-induced de Sitter phase even in 
the case of phenomenologically relevant values of the string coupling, reproducing the standard model couplings today. 
However, in the context of our microscopic model derived from strings, this case would correspond (c/. Eq. ( |68[ )) to a 
large Hubble parameter (and hence curvature) of order M s 2 = 1/a', so our lowest-order approximations would not be 
valid (this of course does not preclude our discussion on the stabilising properties of the condensate from being correct 
in a non-pcrturbativc string context). Hence, it seems that, at least within our string effective context, one cannot 
avoid making the assumption of large negative values of 4>q in the discussion of late-era D-foam vacuum fluctuations 
condensates. 
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2. Matter- dominated phase and statistical contributions of D-particle populations to vector-field condensates 

We next proceed to discuss the massive phase of the recoil vector field, where matter effects dominate, and in this 
sense the contributions of the classical effects of recoil velocity fluctuations, obtained after averaging over D-particle 
populations, to condensates of the vector field strengths dominate over quantum fluctuation effects of D-particles. 

Assuming only electric-field type backgrounds, which is the background considered so far, we do have for such 
classical averages: 



<C F 0v F u > = < E i »= 3 ^ 2 ( f ) > . 
•C F i0 F] > = - < Ei Ej >= -gtj a 2 (t) , 



(78) 



where we assumed isotropic recoil fluctuations for simplicity, and we denoted generically the recoil fluctuations of the 
field strength (averaged over D-particle populations) as & 2 (t), which is a function of the cosmic time. We note that 
scale factor dependent contributions are included, as in Eqs. (49 1, (56 1. The metric gij denotes the FLRW space-time 
background. Using Eqs. (49), (50 1, we get 



a 2 (t)=4H 2 a 2 (t)a 2 =A^-P , H=± 



(79) 



We next notice that, irrespectively of the specific form of the condensate, in case these classical (statistical) effects 
are dominant over the true quantum ones, Eqs. (62), (78) imply 



c*t(i) = -12a 2 (t) < 



a s (t) = +^a t (t) = -4a 2 {t) 



and the corresponding Dirac-Born-Infeld analysis Eq. ( 65 ) yields positive energy density 



DBI-rccoil 



3 A 



d 2 {t) > , 



(80) 



(81) 



but negative pressure 



and an equation of state 



V 



DBI-rccoil ^ ~2/,\ 1 DBI-rccoil 



WBBI- 



(82) 



(83) 



This is the limiting case for inducing acceleration in the universe, which is fine given that we are in the matter- 
dominated phase, where the decelerating effects of matter cannot be ignored. Hence, all this result tells us is that 
the matter-induced recoil effects of D-particles are not sufficient to create a cosmological-constant-type accelerating 
universe, unlike the quantum fluctuation effects 5 . 



However in our model there are many other factors that contribute to the equation of state and to acceleration, bulk D-particles plus 
closed string sector quantum effects contributing to brane tension. Moreover, in our analysis in this work, we have ignored the presence of 
recoil-induced "magnetic-field" contributions. We do notice at this stage that for stability reasons on a D-brane universe, magnetic-field 
type backgrounds are also necessary. In our recoil case this is easily achieved, as discussed in Ref. 1241 . by considering angular momentum 
type contributions in the world-sheet deformations (we give them here in flat space-time for brevity), which arc still compatible with 
the logarithmic conformal algebras, 

^angular momentum ^ [ u ^ t] ^ 0(t - t )d n X k = f f ^ L% k (t - t )8 a X* df,X k + O (S(t - t ) )) , (84) 

and yield magnetic-field type background in target space £?; ~ Uj, with «j the recoil velocity of the defect. Such backgrounds do 
contribute upon performing statistical averages over D-particle populations, terms <K BiBj S> ~ <C EiEj 3> ~ Sija 2 (t), thereby 
making the situation entirely analogous to the one considered in Ref. |23l and described briefly above. The result of these classical 
contributions then in this case is a fluid for the vec tor field describing the dynamics of D-particle-recoil which has positive pressure 
and energy and thus, in the case the constraint Eq. j">2| is ignored, behaves as a relativistic matter fluid with p DBI = (l/3)p DBI > 0. 
We shall not consider, though, such magnetic field contributions in what follows but we felt stating their potentially important role for 
completeness. 
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We will use Eqs. (62), (80) in the equations of motion Eqs. (60) and (59), but keeping the matter stress energy 



tensor terms T™. 



Let us first analyse the graviton equation ( 59 ) . We have 



Rfiu — -j9^v + g^v 87rG c ff A e g — 87rG ff T™ + 87rG c ff 



(85) 



whereby 

G ff 
A cf f 

and 
with (X(x)) 
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Ae* + T 3 e 3 ^ (A(.t)) 
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A Q A" + ^27ra'e-*° 
9s0 



|r 3 |47ra 



7^ [^(1 



^iZ)]'". (86) 



The form of the lagrange multiplier A is found by contracting the vector equation of motion, Eq. ( 58 ) , with A 11 and 
then applying the constraint Eq. ( 52 ) . 



Next, we use the dilaton constraint Eq. (60) to rewrite the term A e 2<6 ° / kq in A e ff. For convenience we split A e g 
into two parts such that 



with A 



A off = AS+A2, 
T 3 e 3 * 



(o) = 

off — 



and A 
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F af3 F r 
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01 



(87) 



This allows the separation of those parts of A e ff which depend on <r 2 terms from those which are cr 2 -independent. 
Similarly, G G ff can also be split as 



and G 



G ff = 




"eff = 


l 


8^ 


^(1) _ 


ae 



.(0) , ^(1) 

"err + u cs 



1 aT^ 
— + — 

K 



9s0 

~32^r 



1 aT^ 

— + — 

K o 9s0 



The constant G^ plays the role of the gravitational constant in this scenario. The background configuration of the 
fields has already been stated in Eq. ( 45 ) and Eq. (50). This then allows us to find the background solutions to the 
equations of motion. 



B. The Background Solution 



For the background, the configuration of the gauge and graviton fields has been explained in the earlier sections 
and was shown to be 



5oo 
9a 

A 

A; 



-1 , 

a 2 (t)Sij , 

-2na/\T 3 \e-+° , , , 



9s0 

-a 2 (t)ui . 



uV(t) 



1/2 



(89) 
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Furthermore, assuming that the matter in our system can be treated as an ideal fluid, it is found that 

T m t = {p + P)w^w„ + P5% , 



(90) 



where p and P are the density and pressure of the fluid and w M is its velocity vector field, normailised to be timelike. 
Thus 



T 







P6 i 



(91) 



With the background configuration defined, we begin looking at the background solutions by examining the vector 
equation Eq. ( 58 1 . By contracting Eq. ( 58 ) with A M is possible to find the form of the Lagrange multiplier in the 
background. This is useful as it can then be used in the other equations of motion when calculating their background 
forms. Thus from Eq. ( 58 ) we can see that 



27ra'|T 3 |e-^ _ 

g s o 2 



[F^(l-ae- 2 ^R)]' v . 



(92) 



With our background and after taking the vev of the recoil velocity in accordance with Eq. (54), it is found that the 
above relation gives 
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(93) 



Next we can see that the vector equation of motion, Eq. (58), is linear in the gauge field. Thus, when looking at the 
spatial component of the equation, all the terms will be proportional to it; and by Eq. ( 54 1 such terms will be zero 
so the spatial component of the equation vanishes. The temporal component remains non zero as Aq from Eq. (89) 
is not proportial to Uj. Thus we find that the temporal component of the vector equation in the background yeilds 

-1 1/2 



= 'dala 2 
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(94) 



Next we can look at the dilaton equation, Eq. ( 60 1 . This gives 
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(95) 



However a more convenient way to use the dilaton equation is to use it as a constraint in the graviton equation, as 
was done in Eq. ( 87 1 . Thus when moving on to the graviton equation we shall be using the dilaton constrained version 
of this equation. When considering the background graviton solution, we can write the modified Einstein equation as 



R 



~ giu ,R = SvrG^ 



[Ml/ 



T: 



DBI 



(96) 



Here Tr^ represents those terms coming from the effects of the Dirac-Born-Infcld field; it is defined to be 



T° BI = 8ttG 



(97) 



Note that for the background solution the (A(x)) term in the definition of A^g is zero upon application of the constraint 
equation. In this case we can write the modified continuity equation as 



T, 



jl V 



T, 



DBI 
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(98) 



For the FLRW background we are considering and denoting by p tot the total density contributions coming from matter 
(p m ), the DBI part (p DBI ) and any dark fluid (p A ), we find, upon assuming a perfect fluid form of stress energy tensor, 
that 



dp tot 
da 
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3P* 



, 



(99) 
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just as in the standard gravitation case. In order to investigate the effects of p DBI in the matter dominated phase of 
structure formation we note the following. The contribution from p A can be neglected in this phase. Moreover, for 
pressureless matter, the total pressure, P tot , has only Dirac-Born-Infeld contributions (including the constraint term 
Eq. (51)) such that P tot = P DBI . The equation of state of the total fluid can then be used leading to 



„DBI 



pt 



(100) 



We thus get the following first order differential equation 



dp b 



da 



with solution 
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where C is a constant to be determined from boundary conditions. The contribution p is 
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(103) 



For the (late) epoch of structure formation, we observe from Eq. (1031 that the terms proportional to CTq ~ a(t) 3 (3 
(cf. Eq. (56)) are suppressed by at least one inverse power of the scale factor a(t) ~ t p . Hence, to leading order, for 



a late era we can write 



DBI 
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(104) 



The background form of the Lagrange multiplier was calculated previously in Eq.(93). Upon examination again we 
see that the Lagrange multiplier terms are suppressed by at least one inverse power of the scale factor. Thus for late 
time eras, the DBI density contribution can be taken to be negligible. 

The total energy density p total = p m + p DBI is determined from Eq. ( 102 ). Matter dominance already required that 
the DBI terms should be suppressed compared to the standard 1/et 3 terms, which we have now explicitly shown is 
the case. Thus, to a good approximation the matter dominated era scale factor, scaling with the FLRW time, can be 
taken to be a ma tter(*) ~ t 2 ^ 3 , as in the standard gravitational case. Lastly then, we can set the constant of integration, 
C, to be equal to the ordinary matter density today, p m Q. We find it convenient to work in units of the critical density 
today p cr = 1, in which p m Q — Q m in the standard notation. Thus, to leading order for the era of structure formation 
where the effect of the cosmological constant can be neglected, 



a, 



(105) 



Some important comments are now in order concerning the nature of the f2 m . In our model, as already mentioned, 
there are contributions to Q m from the D-particles (due to their masses) as well as the conventional dark matter 
components (if there are any). However, due to the competing effects of the nearby and far-away bulk D-particles 
(cf. Eqs (§, gj), there may be significant cancellations in the D-particle contribution to f2 m . As we shall discuss in 
Section |VI| to isolate the effects of D-matter we shall consider the extreme case where f2 m is only baryonic and the 
D-particles play an important role in growth mainly due to their interaction with neutral matter, which provides the 
seed for the appearance of the recoil vector field and its associated perturbations. 

Finally we can write down the explicit form of the density and pressure equations coming from the dilaton con- 
strained graviton equation of motion. Here we shall not assume any particular era of the universe and give the full 
general equations, unlike previously when considering the continuity equation where only the matter domination era 
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was considered and the effect of p\ could be neglected. Below we have also incorporated the form of the Lagrange 
multiplier, Eq. (93). 
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(106) 



Here we can see again that the contribution of those terms proportional to cr 2 are supressed in late eras by at least 
one power of the scale factor. We do notice, howe ver, that these contributions, to leading order in power s of the 



those terms were dominant. The constant terms 



T 3 e 3g ° 
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scale factor, i.e. the terms of the form a^a 2 H 2 in (|106|, do satisfy the equation of state of the DBI fluid (|100|, if 



on the other hand, play the role of the cosmological constant. 



Thus we reproduce, as a further consistency check of our perturbative approach, the standard contributions to the 
density and the pressure in the matter era, obtained previously, with our extra factors coming from the DBI part of 
the action all being suppressed by the scale factor. This allows us to treat their effects as subleading in the late eras 
of the Universe evolution, to which we shall restrict ourselves below. 



C. The Perturbed Solution 



Let us now consider the perturbed metric; it reads 

.goo = -1 + 2*, g i} =a 2 (t)(l + 2<Z>)S l 



The perturbed vector held assumes the form 



A„ 



i{t) (-a(t)tip + A^j , 
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where we used the background solution Eq. ( 50 1 . 



Caution should be exercised when we parametrise the vector held perturbations, which should respect the timedike 



constraint for the vector held Eq. ( 52 1 . One should use the gauge covariance associated with the general coordinate 



diffeomorphisms in order to parametrise the temporal components of the vector perturbations in terms of those of the 
metric [25] . Below we summarise briefly the situation and state the hnal result, relevant to our purposes here. For 
details we refer the reader to Ref. [35]. The reader should recall that, in contrast to the TeVeS models, in our case 
the dilaton is constant and its perturbations are ignored as sub-leading for reasons stated previously. Moreover, the 
spatial components of the background vector field are non vanishing, except when considering averages of Eq. ( 55 ) over 



populations of D-particles in the foam. Nevertheless, these do not affect the general form of the vector perturbations, 
which are of the same form as for TeVeS models. 

The perturbation of the timedike vector field A^, Eq. (1081, in a conformal metric background such as 



.900 = a 2 (v) , 9ij = a 2 {v) 

(where conformal time r\ is related to cosmological time t in the standard way) can be decomposed as follows: 
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where overdine denotes quantities pertinent to the conformal background, is a timedike unit Killing vector tangent 
to the geodesies, t^t^ — — 1, such that = 0, and q^ v = 5^ + t^t u is a projector appropriate for the conformal 
metric [3S], such that Tt v A v = 0, 5'a? < t = Q^v The fields £ and S are scalar modes, and /3 M are vector modes such that 
lifiu = 0. There are two independent vector modes, as a result of the constraint Eq. (52 1 and gauge invariance 
under diffeomorphisms. The scalar mode C satisfies [25] 
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There is gauge freedom in transforming the various fields under general coordinate transformations, which results 
in simplifications (upon gauge fixing) of the vector perturbations. To this end, let one consider an infinitesimal 
diffeomorphism £ M = where a is the scale factor of the conformal metric. The vector £^ can be decomposed as 
£ M = — £t M + g l/ M V„V + ' i V> with f 1 ^^ = 0, which contains two scalar modes £ = and ip, such that \/ 2 ip = ^"V,,^, 
and two vector modes with the property q^V = 0. Under such diffeomorphisms, the vector perturbation 
transforms as 



A-n — 



1 



(112) 



we recall that we ignore dilaton perturbations, which would in general contribute. 
From Eqs. (110) and (112) we observe that the scalar mode £ transforms as 

1 



the scalar s transforms as 



c = c + -e , 

a 



a 



(113) 
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while the vector modes /? M remain gauge invariant. 

Using the gauge freedom endowed in Eqs. ( 113 ), ( 114 ) we may then gauge fix the perturbations of the vector field, 
by making the gauge choice [25] 
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where \& is defined in the metric perturbation Eq. (107). Thus, the form of the (gauge fixed) vector perturbation 
reads 



with A satisfying Eq. (Ill), or cquivalently 
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setting constants to zero. A shift back to the Robertson- Walker metric using a standard co-ordinate transformation 
leads to 
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Here the 
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factor is included from dimensional considerations. The last perturbations which need to 



be defined are the perturbations of the matter energy momentum tensor. In our system Eq.(90) can be perturbed to 
find that 
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We can now begin to examine the perturbed equations of motion, beginning with the graviton equation. In doing so, 
it is useful to move into Fourier space, such that we get the following conversion 

di = iki . (120) 

This converts the partial differential equations into ordinary ones, making them easier to manipulate. The 00 com- 
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ponent of the graviton equation yields the following 
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where k 2 = kik 1 is the amplitude of the spatial components of the momentum scale. Similarly, the ij component of 
the graviton equation yields 



1 

,2„2 



9s0 



+2er^a 2 $o\ < H 



+2cr 2 a 2 *(5 



3 a 



6ae 



5H 2 



a* 



2H 2 



~20o 



2a 2 



+ 



4a 4 



a^ a 
+ 6ae" 2 ^ 



2iTa 6i? 2 a 
+ 



H 4 



4a 
a 

7H 2 k 2 



W + 2i? 3$ + * 



2$ 



a 

a 2 ZTa 

^2 ^ 



6a 2 



a 



12cr 2 a 2 ae- z<?>0 $<5 1 , <j 4# 



20o ( T ) \i 



Fa 
a 



2ff 4 4 

12CT 2 a 2 ae~ 2 ^*^ jsff 3 



fc 2 


3a, 


12a 4 


a 


5Hd\ 


a 





H 2 (3 + 10a 2 



-12cr 2 a 2 ae- 200 ($ - *J ^ + 6cr 2 i^ae -2 * - *) 

fc 2 <5\- 



2a 



1 _ 6ae" 2 ^ ( if 2 + 

a 



(122) 



In the above we have separated out those components which come from the standard Einstein Hilbert part of the 
action from those terms which appear due to the effect of the D-particles. The standard terms are those which are 
proportional to 1/kq + aT^e^ /g S Q. At this point a very useful simplification can be made. The T m% j equation has 
two types of terms, those proportional to 5 % j and those proportional to k % kj. By contracting the equation with the 
following tensor operator 



1 



S J ., 



(123) 



all the terms proportional to S l j, including those on the left hand side of the equation are set to 0. Thus, after this 
operation on Eq. ( 122 1 we find that 



* - $ = 



(124) 



Thus, in our system, just as in the standard gravitational case, $ = ^ is a valid solution which we shall now be using 
throughout. As we shall show later on, the perturbations of the recoil vector field do exhibit growing modes and 
thus can participate in late-era structure formation, unlike the case of TeVeS-like models in [TT], where there is a 
non-trivial difference ^ — $ ^ that is considered as the main source of the growing mode [26] . 
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Upon using Eqs. ( 124 f , (1191 and (1051, we obtain for the dimensionless over-density parameter 5p/p the result: 
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From this equation we can calculate the growth of density perturbations, provided that the evolution of the metric 
and vector perturbations are known. An important point to note is that the vector perturbation terms, £ and £ only 
appear due to the {\{x))A fl A l , term of T^v in Eq. (87). Thus the constraint term is vital in coupling the vector 
perturbation to the density perturbations. 



VI. GROWING MODES AND LARGE-STRUCTURE FORMATION DUE TO D- PARTICLES 



We now proceed to calculate the evolution of the vector field perturbations. We notice that the remaining equations 
of motion provide the constraints which allow the evolution of these perturbations to be calculated. Under the 
perturbations of the metric it is noticed that, since all the terms in the vector equation are linear in the gauge field, 
any perturbations which appear from the metric will always come with a Uj term. Such terms, when averaged over, 
are zero (as discussed previously) . Th us, in our case metric perturbations only appear in the vector equation through 
the presence of W = $ in A (cf. Eq. (118)). The perturbed vector equation reads 
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where 
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where we remind the reader that (X(x)) is the background Lagrange multiplier, the form of which is given in Eq. (93 1 . 
In Eq. ( 126 ) the signs of terms &i and 6 2 play a crucial role in determining if the vector modes of the perturbation 
enter a growing mode or a decaying one, a result which is analogous to that found in Ref. |llj . In general, growth 
will be seen when the contribution from 62 is negative. 

To this end, we first note that for large epochs of the Universe, namely the radiation and matter dominated eras, 
it can be readily seen that (X(x)) > 0. In our analysis below we shall assume values of <fio = 0(1), which may be 
phenomenologically desirable, since the string coupling g s — determines the gauge couplings of the low-energy 
theory. Looking at &2 it can be seen that in this case the (A(x)) term in crucial in determining if 62 will change sign 6 . 



6 We note for completeness that in models with large and negative values of (j>g, the terms — ^M- 2.1 in fe 2 , although sub- leading 

at late eras of the Universe, since they contain higher orders of the scale factor and its derivatives, nevertheless could have non-negligible 
contributions to 62 due to the denominator, which may be sufficiently small for sufficiently large and negative (f>o, such that the over all 
contribution of this term to 62 is significant. In our case with 4>0 = 0(1), though, such terms will not play any role in our analysis. 
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FIG. 2: Vector perturbation ( as a function of the scale factor a(t) in the matter-dominated era, for different values of the 
D-particle recoil velocity variance a\ oc /3/a 3 . The string mass M s is assumed to be 10 TeV, c/>o = 1 and k = 0.5 Mpc -1 . 



Since (A(x)) is proportional to the variance of the recoil velocity (cf. (93)), <7q, we see that it is the magnitude of c 2 , 



which will determine if the vector perturbations will enter a growing mode. 
Finally we examine the perturbed dilaton equation. We get 
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We remind the reader that in our approach we assumed that any dilaton couplings to matter have been suppressed, 
and this is the reason why no matter contributions appear in the perturbations associated with the dilaton constraint 
Eq. ( 128 ). Here we can see that the dilaton equation specifies the evolution of the metric perturbations entirely. With 
the form of $ then inserted in to Eq. (126), the form of £ is found, and finally with both ( and <3? used in Eq. ( |125[ ), 



dp/p can be calculated. 

The numerical results for the matter dominated era, that is a scale factor a ~ t ' , are shown in Figs. §and[3j 
which also shows the effect of altering the D-particle recoil velocity variance, (cf. Eq. (56)). A baryon only scenario 



was considered for the numerical results as an extreme case, in order to isolate the recoil-velocity effects of D-matter. 
From the figures it can be seen that the magnitude of the variance of the recoil velocity of the D-particles plays 
the crucial role of allowing matter density perturbations to grow sufficiently to allow for structure formation. In 
particular, Fig. [3] shows that when the parameter /3 is sufficiently small, j3 < O(10~ 4 ), then Sp/p begins to show 
oscillations and no longer exhibits growth. This is to be expected as in the standard gravitational paradigm these 
same oscillations appear in the absence of dark matter. The role of dark matter is in our case being played by the 
coupling of the vector field perturbations to Sp/p. Thus for higher values of /3 the vector perturbations are strong 
enough that they can drive structure formation in much the same way that dark matter does in the standard case. 
Our point is not, however, to do away with dark matter especially as dark matter candidates fall out naturally from 
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FIG. 3: Over-density parameter |<5| = 5p/p as a function of the scale factor a(t) in the matter-dominated era, for different values 
of the D-particle recoil velocity variance a\ oc /9/a 3 . The string mass M s is assumed to be 10 TeV, 4>o = 1 and k = 0.5 Mpc -1 . 
\5\ is taken to be 10 -5 today. 



the string theory framework our approach relies upon. The point is rather that there may be alternative methods of 
driving structure formation which are in operation in conjunction with dark matter; the vector perturbations of the 
approach presented here are one such candidate which has now been shown to have the required characteristics to 
drive structure formation. 

We can also note that the theory is largely insensitive to the string mass, with the critical value of (3 to induce the 
growing mode remaining approximately of the same order over the range 10 TeV < M s < 10 18 GeV. This can be 
readily seen by looking once more at the expression for the coefficient 6 2 appearing in the perturbed vector equation, 



Eq. (126). As we noted earlier, only when <fio is large and negative will those terms proportional to a = 1/Mf play 
a significant role. In the above numerical analysis we took <fio = 1, thus these terms have negligible effect. In the 
background form of the Lagrange multiplier (A(cc)), given in Eq. ([93]), we see that the string mass dependence there 



is actually cancelled by the string mass term appearing in the definition of a 2 given in Eq. (56). Thus (A(a;)) is 
independent of the string mass and so the value of (3 for which the vector field will enter its growing mode is also 
independent of the string mass. 

A further point of interest is the behaviour of the perturbed vector equation in the radiation dominated era. In the 
radiation dominated era the scale factor a ~ t 1 ! 2 and as result, 



H 2 + - = 
a 

2H 3 ~ — -'^=0 (129) 



a a 



This can be used in the expression for the coefficient &2 in Eq.(127). 62 now becomes, 



62 = -2(A(ar)) . (130) 
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Thus in the radiation dominated era &2 is always negative, since (\{x)) > (cf. Eq. (93)), which as mentioned before 
is the condition for the growth of the vector perturbations. Thus in this era the vector perturbations are always 
in a growing mode, irrespective of the specific density of the defects, provided of course the latter is non zero and 
sufficiently large so that the above formalism based on the existence of a recoil-vector field is valid. Once again then 
the vector perturbations show a behaviour similar to that of dark matter by being able to grow during the radiation 
dominated era when baryon density perturbation would not be able to cluster due to the radiation pressure. This 
allows the seeding of the density perturbations which would drive structure formation to begin early in the history of 
the universe. As the brane Universe enters the matter-dominated era, the growing mode persists only above a critical 
density of defects, as we have discussed previously (cf. Figs. 2]3), which is an interesting feature of the model. Of 



course in the context of our string models, in this era, other (conventional from a particle physics viewpoint) candidates 
for dark matter, such as supersymmetric partners of mater excitations or gravitinos do exist and contribute to the 
dark matter spectra. However, in view of the assumption of sufficiently dense populations of the defects on the brane, 
required for the "medium" interpretation, the cosmological constraints of such models, especially as far as collider 
searches for supersymmetry are concerned [27], are expected to be modified, depending on the mass range of the 
D-particles [7]. 

VII. CONCLUSIONS AND OUTLOOK 

In this work, we have analysed some microscopic string theory models of modified gravity arising in the low-energy 
limit of brane worlds containing space-time point-like brane defects (D-particles). Dense media of such D-particles 
can exist at early eras of the Universe, without the danger of over-closing the Universe, for specifically stringy 
reasons. Propagation of neutral matter, such as neutrinos, on such backgrounds leads to effective gravitational theories 
containing, in addition to the traditional graviton and dilaton fields of the gravitational multiplet of the string, also 
vector gauge fields, describing the recoil of the brane defects during their topologically non-trivial scattering with the 
string matter. The vector fields are associated with the recoil velocities of the defects and as such satisfy a given 



constraint Eq. (52). The presence of a recoil velocity locally breaks Lorentz invariance, which however is assumed 
to be restored on average over large populations of D-particles, since the velocity expectation value vanishes, leaving 
only the variances of the velocities to be non zero. 

The effective Lagrangian describing the low-energy dynamics of this model contains a Dirac-Born-Infeld type of 
lagrangian for the vector field, coupled non-trivially to space-time curvature. By solving the associated equations of 
motion we have determined our background configuration, over which perturbations were considered. We have also 
demonstrated the consistency of the solutions with the conformal invariance conditions of the associated stringy a- 
model. The equation of state of the Dirac-Born-Infeld fluid of the vector field has also been considered, and this result 
can be derived exactly to all orders in a 1 . The gravitational sector of the model, unlike the vector field, cannot be 
studied exactly but only perturbatively in powers of a' . In this work we restricted ourselves to considering space-time 
curvature terms in the effective action up to order O (a') , which suffices for our low-energy considerations at late 
epochs of the Universe. 

By considering perturbations of the vector field, we have shown the possibility of a growing mode in the matter 
era, for sufficiently large values of the variance of the recoil velocities of the D-particle. The mode also exists in 



the radiation era. The constraint Eq. ( 52 ) of the vector fields is essential in ensuring that the growing mode exists. 
This allows the seeding of the density perturbations which would drive large-structure formation to begin early in 
the history of the universe. This feature is shared by TeVeS models in alternative to Dark matter scenarios, but 
the main difference of our model lies on the fact that in our case the growing mode is independent of the difference 
of the two gravitational potentials \& — <f> which vanishes here. This is consistent with the equations of motion for 
the graviton and dilaton fields, and thus the associated conformal invariance conditions of the stringy er-model. In 
addition, although in our theory, D-matter plays a role analogous to dark matter as far as large-structure formation 
is concerned, this is only one component, given that the underlying superstring inspired models do involve additional 
components coming from the supersymmetric partners of the standard model particle sector of the theory. In this 
sense, the phenomenology of these models is different from conventional low-energy string-inspired effective theories 
and can depend crucially on the density of D-particle defects (and the magnitude of their masses M s /g s ) in the current 
era. In this spirit, comparison of our class of models against gravitational lensing data is essential in determining 
the amount of dark matter present in the centre of galaxies, where the concentration of D-matter is expected to be 
significant. 

In our solutions in this article we assumed that the dilaton field was constant. Extensions of our model to incorporate 
time-dependent dilaton fields are envisaged, although in such a case the reconciliation of the model with the particle 
physics phenomenology may be subtle, given that the exponential of the dilaton is connected to the string coupling, 
and the latter to gauge couplings of the low energy field theories coming from the string. If we want standard model 
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physics to be reproduced today and the findings of the theory of Big Bang Nucleosynthesis not to be disturbed, then 
one should apply stringent constraints on the varying dilaton fields. Moreover, non constant dilaton fields may have 
significant impact on the growth of galaxies, and thus severe constraints from the relevant astrophysical data are likely 
to be imposed, which can easily rule out most of such cosmologies [2"8] , 

Finally, we would like to close by stressing the fact that the detailed cosmology of such models is still in its infancy, 
but we believe that the (rather toy) model we studied here exhibited several interesting and non-trivial features 
that are worthy of further detailed investigations in more realistic models and in very early eras of the Universe. 
An interesting aspect we would like to study is the role of the D-particles and the associated vector recoil field in 
inflation. Such issues may be tackled by considering colliding branes in our context (c/. Fig. [lj. Brane collision would 
play the role of a cosmically catastrophic event. During the collision, the concentration of massive bulk D-particles 
that are trapped in the bulk space between the colliding brane worlds can increase in such a way so that they can 
collapse to form black holes. The latter would then evaporate, with the standard model particles being trapped on 
the brane. The gravitational excitations of the string, on the other hand, are capable of escaping in the bulk regions. 
The collision process can lead to inflationary expansion on the brane universe [29], while the black hole evaporation 
can contribute to reheating of the inflationary universe and subsequent graceful exit from the inflationary phase. The 
particular role of the recoil vector field during the interaction of string matter with the defects can be investigated 
at such very early epochs of the Universe. Although at present these considerations appear to be mere speculations, 
nevertheless we believe that these issues can be tackled at some detail within the framework of our string-modified 
gravity. We hope to be able to report progress on some of these topics in a future publication. 

A final comment we wish to make concerns the density of D-particles on the brane world. Although cosmologically 
the later is left unconstrained, as we have already mentioned, given that the issue of overclosure of the brane Universe 
does not arise, nevertheless, the presence of a finite density medium of D-particles on the brane, with which neutral 
particles such as neutrinos and photons interact non trivially, implies non-trivial "optical" properties for the brane, 
such as a refractive index. The latter will manifest itself as energy-dependent delays in the arrival times of photons 
emitted simultaneously from a high-energy astrophysical source, such as a Gamma Ray Burst or Active Galactic 
Nucleus. As discussed in |14[ [30] . during each interaction of a high energy neutral particle (represented by a string) 
with a D-particle, there is a delay in the re-emission process of the particle due to capture which is linear in the 
incident particle energy E, St ~ E/Mg, where M s is the string scale (we work in units of c = 1). The total delay turns 
out to be proportional to the number of D-particle defects per string length, rj£ 3 , encountered by the propagating 
photon, times the distance travelled, Ai tota i ~ j^EL. This implies an effective suppression scale M eS = M s /?y 4 . 
The sensitivity of current observations of photons from sources at redshifts z < 3 [141 130] to such linearly suppressed 
delays is such that M is of order of the Planck scale (10 19 GeV). Since from particle physics experiments we know 
that M s > O(10) TeV, we observe that an 7/4 as low as 10 -15 can lead to observable effects from refractive index 
measurements for low string scales. In our models discussed here, we may easily imagine a depletion of defects in 
the bulk during epochs z < 5, so that their densities on the branes fall below such values, thereby avoiding any 
constraint from the above refractive index measurements. If our defects on the brane (D-matter) are one or two 
orders of magnitude more dense that the standard (conventional) dark matter, which is a quite natural situation, 
these optical constraints will still be negligible. Finally, we remark that the presence of D-matter may affect the 
peaks of the Cosmic Microwave Background Radiation, depending on the respective densities at that epoch, but this 
analysis requires a separate study, which goes beyond the purpose of our present article. 
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